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In an extension of the Nambu-Jona-Lasinio model where the quarks interact with the temporal 
giuon field, represented by the Polyakov loop, we explore the relation between the deconfinement and 
chiral phase transitions. The effect of Polyakov loop dynamics on thermodynamic quantities, on the 
phase structure at finite temperature and baryon density and on various susceptibilities is presented. 
Particular emphasis is put on the behavior and properties of the fluctuations of the (approximate) 
order parameters and their dependence on temperature and net-quark number density. We also 
discuss how the phase structure of the model is influenced by the coupling of the quarks to the 
Polyakov loop. 



PACS numbers: 25.75.Nq, 12.39.Fe 



1. INTRODUCTION 

Quantum Chromodynamics (QCD) exhibits dynami- 
cal chiral symmetry breaking and confinement. Both 
features are related with global symmetries of the QCD 
Lagrangian. However, the relation of spontaneous chi- 
ral symmetry breaking and confinement still remains 
an open issue. While the chiral symmetry is exact in 
the limit of massless quarks the Z{3) center symme- 
try, which governs confinement, is exact in the oppo- 
site limit, i.e., for infinitely heavy quarks. In order to 
obtain a unified picture of confinement and chiral sym- 
metry breaking several effective chiral models have been 
studied P, i, [1 i, i, H, 0, H, Hi ■ Recently an extension of 
the Nambu-Jona-Lasinio (NJL) model [lO| . | T l| was pro- 
posed and developed 0, S S [13, Il4l. llSj to address 
this question. 

The NJL model describes interactions of constituent 
quark fields. It exhibits a global SU(3) symmetry that is 
a replacement of a local gauge SU(3) color transforma- 
tion of the QCD Lagrangian. Thus, color confinement is 
lost in the NJL dynamics. Recently, color degrees of free- 
dom were introduced in the NJL Lagrangian through an 
effective gluon potential expressed in terms of Polyakov 
loops Such a potential was constructed to pre- 

serve the Z{3) symmetry of the gluon part of the QCD 
Lagrangian and it has its origin in the recently devel- 
oped effective models with Polyakov loops as dynamical 
fields [IB]- An extended NJL Lagrangian (PNJL model) 
contains unified properties of QCD related with Z{3) and 
the chiral symmetries. The interactions of quarks with 
the effective gluon fields in the PNJL model is included 
through covariant derivatives. Furthermore, due to the 
symmetries of the Lagrangian, the PNJL model belongs 
to the same universality class as that expected for QCD. 
Thus, such a model can be considered as a testing ground 
for studying the phase structure and critical phenomena 
related with the deconfinement and chiral phase tran- 
sitions. This is particularly interesting since there are 
still limitations in the applicability of lattice gauge the- 
ory (LGT) to QCD thermodynamics at large net quark 



densities. 

Recently, it was shown [1, [H, [H, ^ that the PNJL 
model, formulated at finite temperature and finite quark 
chemical potential, reproduces some of the thermody- 
namical observables computed within LGT. The prop- 
erties of the equation of state Q, the in-medium modi- 
fication of meson masses [l^ as well as the validity and 
applicability of the Taylor expansion in quark chemical 
potential used in LGT were recently addressed within the 
PNJL model [l3,[ll,[Il- In Ref. [ll the model was ex- 
tended to a system with finite isospin chemical potential 
and pion condensation was studied. 

In this paper we will consider the fluctuations in var- 
ious channels at points in the phase diagram near the 
phase boundary. In particular, we focus on the behavior 
of the net -quark number fluctuations and the suscepti- 
bilities of the order parameters'*^^. The susceptibilities of 
the Polyakov loop and its conjugate as well as the chiral 
condensate will be introduced and analyzed. The rela- 
tion between the chiral and deconfinement phase transi- 
tions will be quantified through the susceptibilities. Our 
calculations are performed within the mean field approx- 
imation and in the PNJL model, including a non-local 
four-fermion interactions to regulate the divergent mo- 
mentum integrals. 

The paper is organized as follows: In Section [5] the 
PNJL model is formulated and the relevant thermody- 
namic potentials are derived. In Section [3] the influence 
of the Polyakov loop dynamics on various thermodynam- 
ical quantities is discussed. In Section 2] and Section [S] 
we introduce the susceptibilities of the order parameters 



#^ Wc will be somewhat lax, and refer to the Polyakov loop and 
the quark condensate as order parameters, although, due to the 
explicit breaking of the symmetries, neither of them is strictly 
speaking an order parameter. However, as found in LGT calcu- 
lations, both are useful quantities for identifying the location of 
the phase boundary. Thus, it seems that the explicit breaking 
of both the chiral and center symmetry in QCD is in this sense 
weak. 
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and focus on their properties. Concluding remarks and 
discussion are presented in Section [51 



2. TWO-FLAVOR NJL MODEL WITH THE 
POLYAKOV LOOP 

Various methods have been proposed to account for 
interactions with the color gauge field in effective chiral 
models P, d, d, H, 0, 01 ■ One of these is the extension of 
the Nambu-Jona-Lasinio (NJL) Lagrangian by coupling 
the quarks to a uniform temporal background gauge field, 
which manifests itself entirely in the Polyakov loop 0, 
H, [1] . The PNJL Lagrangian for three colors {Nc = 3) 
and two flavors {Nf ~ 2) with non-local four-fermion 
interactions is given by 

C = 4>{ip-m)i: + '4!iljo^/j-UmA],^A];T) 
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iqix)qix)r + {q{x)l-i^fq{x)Y , (2.1) 



where rh = diag(m„, m^) is the current quark mass, jl = 
diag(/i„, ^d) is the quark chemical potential and r are the 
Pauli matrices. We assume isospin symmetry and take 
niu = md = mo and fiu ^ fJ-d ^ A*- 

The Lagrangian is formulated with non-local interac- 
tions, which are controlled by a form factor. This feature 
of the model is implemented in order to deal with the 
ultraviolet singularities that appear in the loop integra- 
tions. In coordinate space, the form factor F{x) for the 
non-local current-current interaction reads: 



q{x) = / d'^yF{x - y)i}{y) . 



(2.2) 

A possible choice for the regulator is in momentum space 
given by p^ : 



1 



1 + {p/kf 



(2.3) 



where f{p) is the Fourier transform of the form factor 
F[x) and p is the three-momentum. 

The strength of the interaction among constituent 
quarks in (|2.ip is parameterized by the coupling constant 
Gs which carries the dimension of a length squared. In 
the pure NJL sector, the model is controlled by four pa- 
rameters: the coupling constant G5, the current quark 
mass TTiQ and the constants a and A, which characterize 
the range of the non-locality. These parameters are de- 
termined in vacuum, for a given a, by requiring that the 
experimental values of the pion decay constant ~ 92.4 
MeV and the pion mass m.^ ~ 135 MeV as well as the 
dynamical quark mass Mp=o = 335 MeV are reproduced. 
For a = 10 the model parameters are A = 684.2 MeV, 
GsA^ = 4.66, Too = 4.46 MeV. The corresponding value 
of the quark condensate is {■tptpy^^ = —256.2 MeV fl^ . 
The relevant parameters of the model used in our calcu- 
lations are summarized in Table HI 



The interaction between the effective gluon field and 
the quarks is in the PNJL Lagrangian implemented (|2.ip 
by means of a covariant derivative 



A^, = 6^oA° , 



(2.4) 



where we introduce the standard short-hand notation 
Af^ = gA^j^^. Here g is the color SU(3) gauge coupHng 
constant and A° are the Gell-Mann matrices. 

The effective potential U of the gluon field in (|2.1[1 is 
expressed in terms of the traced Polyakov loop $ and its 
conjugate $ 



— Tr.Lt 

Nr. 



(2.5) 



where i is a matrix in color space related to the gauge 
field by 



L{x) = V exp 



I I drA^ (x, r) 




(2.6) 



with V being the path (Euclidean time) ordering, and 
P = l/T with A4 = iAo. 

In the heavy quark mass limit QCD has the Z(3) cen- 
ter symmetry which is spontaneously broken in the high- 
temperature phase. The thermal expectation value of the 
Polyakov loop ($) acts as an order parameter of the Z(3) 
symmetry. Consequently, ($) = at low temperatures 
in the confined phase and ($) ^ at high temperatures 
corresponding to the deconfined phase. For the SUc(3) 
color gauge group the Polyakov loop matrix L satisfies 
LL^ = 1, detL = 1 and can be written in diagonal form 



L = diag(e"^,e''^ 



(2.7) 



In general $ in Eq. (12. 5|) is not identical to $. 

The effective potential W($, $) of the gluon field is 
expressed in terms of the Polyakov loops so as to preserve 



the Z{S) symmetry of the pure gauge theory 16 
adopt an effective potential of the following form 



We 
#2. 



z^($,l>;r) _ &2(r) 



rp4 
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(2., 



with 



b2iT) = ao+a, ('^] +aJ^y +aJ^y . (2.9) 



The presence of quarks will clearly affect the Polyakov loop dy- 
namics. Thus, one may ask why this is not reflected in the effec- 
tive potential U. However, at least partly such effects arc treated 
explicitly in the PNJL model through the coupling to the quarks. 
Thus, in order to minimize the double counting problem, we keep 
the Z{3) symmetric form of the effective potential, but consider 
the possibility that the parameters, in particular Tq, may de- 
pend on the number of dynamical flavors Nf, thus accounting 
for quark loop effects not included in the mean-fleld approxima- 
tion. 
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The coefficients Tq, and bi are fixed by requiring that 
the equation of state obtained in pure gauge theory on 
the lattice is reproduced. In particular, at Tq = 270 MeV 
the model reproduces the first order deconfinement phase 
transition of the pure gauge theory. The parameters are 
listed in Table HD 

In the mean field approximation the Lagrangian (|2.ip 
is rewritten as 



where 



Gs 



((g(a;)g(x)))^-Z^(<f>,$;r), 



(2.10) 



where 



yi) 



X 

+ Fix - yi)F{x - y2){q{x)q{x)) . (2.11) 

In momentum space the dynamical quark mass is deter- 
mined by the current quark mass mg, the quark momen- 
tum distribution function f{p) and the chiral condensate 

im) 

Mp=^mo + {M-mo)f{p), 

M = mo-Gs{qq). (2.12) 

where M denotes the dynamical quark mass at p = 0. In 
a uniform system, the chiral condensate 



qq) = / d^pfip)^mp)i'ip)) 



(2.13) 



is independent of the position. 

From the mean field Lagrangian (|2.10p one obtains the 
thermodynamic potential in the following form 

0(M,$,l>;r,//) =Z^($,$;T)H ("^o - M)^ 



2Gs 



6Nf 



-0^ [E, ~ EpiM, = mo)] 



2NfT 
+ Ti-r In 



In 



1 + Le 



-E'-+yT 



-E^-'i IT 



(2.14) 

where we have introduced E'^^^ = ^ /i for the par- 
ticle (-1-) and anti-particle (— ) with Ey = y^jpl^ + A/^ 

being a quasiparticlc energy. The third term in the ther- 
modynamic potential (|2.14p . which corresponds to the 
vacuum contribution, is renormalized by subtracting the 
term Ep{Mp = mo) under the momentum integral. This 
removes the divergence from and introduces an irrelevant 
constant shift of fl. 

Furthermore, by taking the trace in color space, one 
obtains the final expression for the thermodynamic po- 
tential 

(mo - M)2 



17(M, $, T, /x) = W($,$;T) 
d^p 



6Nf 



2NfT 



27r)3 
d'^p ( 
(2^1 



2Gs 

[Ej, - Ep{Mp = mo)] 

In c/W(Af,$,$;r,/x;p) 



In 



(M,$,$;T,Ai;p)]} 



,g(+)(A/,$,$;T,/.;p) 



7(-)(A/,$,|.;r,M;p) 



1 + 3 $ + $e" 



E'--^ /t\ „-E'--^ it I -■iE'--^ IT 



(2.16) 



Under the transformation /i —p, the role of quarks 
and antiquarks is exchanged. Inspection of (|2.16p shows 
that in the PNJL model the charge conjugation trans- 
formation also exchanges the role of the Polyakov loop 
and its conjugate. This is reflected in the relation 
<?(+) (M, $, $; T, p) = g(-) (Af, $, $; T, p) . 

Although the Polyakov loop is complex, the thermo- 
dynamic potential is real. Due to the symmetry in color 
space, the imaginary part of f2 vanishes after performing 
the functional integral [see Appendix [A] . 

An interesting feature of the PNJL model described 
by the thermodynamic potential (j2.15|) is the qualita- 
tive behavior in the low temperature phase. In the 
limit of <I>, $ 0, which is expected at low temper- 
atures, the contribution of one- and two-quark states 
to g^^^ are suppressed and only the three-quark term 
- exp{-3E<'^^ /T) survives. In this sense the PNJL 
model mimics the confinement of quarks within three- 
quark states. The suppression of quark degrees of free- 
dom at low temperatures is, on a qualitative level, sim- 
ilar to confinement in QCD thermodynamics. Thus, 
the PNJL model is better suited for describing the low 
temperature QCD phase than the standard NJL model, 
where the constituent quarks are abundant also at low 
temperatures. However, at least in the mean-field ap- 
proximation, the model only has the three-quark states, 
but there arc no one- and two-quark states, which also 
play an important role at low temperatures. 

In the mean field approximation the dynamical quark 
mass M and the expectation values of the Polyakov loop 
$ and $ are obtained from the stationarity conditions 



(2.17) 



obtained by extremizing the thermodynamic potential 
with respect to M, $ and $. 



(2.15) 
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These conditions yield the foUowing set of coupled gap equations: 



M = mo + 6 GsNf 



(2^)3 Ep 



62(r)$ + 63$^-64($$)$ = 
62(r)$ + 63$^ -64$($$) = 



.g(+)(M,<i>,$;T,Ai;p) 



127V/ f (Pp 



T3 y (27r)3 
l2Nf f (fp 



5(-)(Af,$,$;r,/i;p) 



g(+) (M, $, $; T, p) g(-) (A/, $, T, p) 



T3 7 (27r)3 



,-2£;<+Vr 



-E^-^ IT 



g(+) (M, $, $; T, p) g(-) (A/, $, $; T, p) 



(2.18) 
(2.19) 
(2.20) 



As noted above, the role of <& and $ are exchanged un- 
der charge conjugation (/i — > — /i). This symmetry is re- 
flected in the relation between the gap equations (|2.19p 
and (|2.20p under the charge conjugation transformation. 
In particular, this relation implies that $ = $ for /i = 0. 
Furthermore, we note that for $ = $ = 1, Eq. (|2.18p is 
reduced to the gap equation of the standard NJL model 
without any coupling to the color SUc(3) gauge field. 



3. THERMODYNAMIC QUANTITIES AT 
FINITE QUARK CHEMICAL POTENTIALS 

The thermodynamics of the PNJL model in the mean- 
field approximation is characterized by the potential 
fi($,$,M) introduced in Eq. ((2l^ . In Figs, [l] and [2] 
we show the effective potential of the Polyakov loop U as 
well as the PNJL thermodynamic potential fl in the chi- 
ral limit, at a high temperatures where chiral symmetry 
is restored and the dynamical quark mass M vanishes. 
Consequently, only the Polyakov loop and its conjugate 
are the relevant classical fields in the problem. 

For vanishing quark chemical potential, $ and are 
equivalent. Hence, the thermodynamic potential is char- 
acterized by only one variable that for a given temper- 
ature is determined by the gap equation (|2.19p . In Fig. 
[1] we show the $ dependence of U. This potential ex- 
hibits the expected structure in the Z{3) symmetry bro- 
ken phase with a minimum at a finite value of $ below 
unity. 

The infiuence of quarks on the PNJL model thermody- 
namics is illustrated in Fig. [21 where the ^-dependence 
of n is shown for the same temperature, T = 0.5 GeV. A 
comparison of Figs. [T] and [5] clearly shows that the inter- 
actions of the effective gluon field with the quarks leads 
to a shift of the minimum to larger values of <&. At high 
temperatures the minimum corresponds to $ > 1. 

The dynamical quark mass M and the traced Polyakov 
loop $ are the order parameters of the chiral and Z{3) 
symmetries respectively. Thus, the T and ^ dependence 
of these order parameters is used to identify the phase 
boundaries of the model. In Figs. [3] and 2] we show <!> and 



U/T* 




FIG. 1: The Polyakov loop effective potential U/T^ as a func- 
tion of $ at T = 500 MeV in the chiral limit. 




FIG. 2: The thermodynamic potential Sl/T"* as a function of 
$ at T = 500 MeV in the chiral limit. 



M ai II = 0, as functions of T in the chiral limit. The 
phase change of the model is clearly indicated by rapid 
changes of the order parameters. The Polyakov loop po- 
tential U introduced in the PNJL Lagrangian preserves 
the invariance under the Z{i) symmetry. However, due 
to the interactions with the quarks this symmetry is ex- 
plicitly broken in the model. Thus, the Polyakov loop is 
not an order parameter in the strict sense and the tran- 
sition is a rapid cross over. This is a remnant of the 
" deconfinement" phase transition of the pure Polyakov 
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FIG. 3: Expectation value of the traced Polyakov loop $ in 
the chiral limit as a function of temperature T for vanishing 
quark chemical potential. 
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FIG. 4: The dynamical quark mass M in the chiral limit as 
a function of temperature T for a vanishing quark chemical 
potential. The dash-dotted line denotes the result obtained 
in the NJL model. 



loop model. On the other hand, in the chiral limit the 
chiral transition is a true second order phase transition, 
with an order parameter M, which is strictly zero for 
temperatures above T 242 MeV. 

In Fig. m the dynamical quark mass M obtained in 
the PNJL model is also compared with that of the NJL 
model . It is clear that the coupling of the quarks to 
the effective gluon field shifts the chiral phase transition 
temperature to the higher value. This is a consequence 
of attractive interactions; in the presence of the Polyakov 
loop there is stronger " binding" of the constituent quarks 
in the chirally broken phase. 

The non-locality of the quark interactions introduced 
by the form factor F{x) also affects the critical tempera- 
ture of the chiral phase transition. Suppressing the form 
factor in the PNJL Lagrangian results in a lower tran- 
sition temperature. With the parameters used in our 
actual calculations, the reduction of Tc is on the order of 



In the NJL model, without Polyakov loops, we use the local 
version of the cutoff, i.e., a sharp cutoff. The cutoff is not imple- 
mented in the thermal part, which is anyway convergent. 



15 MeV, if in the thermal part of the PNJL model the 
momentum cutoff is not implemented. 

The introduction of a finite quark chemical potential is 
expected to modify the behavior of the order parameters. 
The n dependence of A/, $ and $ are shown in Figs. [5] 
and [HI There is a clear shift of the chiral transition to 
lower temperatures with increasing ^ as is seen in Fig. [SI 
At a non-zero quark chemical potential the charge con- 
jugation symmetry is broken, which leads to a splitting 
between 4> and As seen in Fig. the Polyakov loop 
$ is decreasing whereas $ is increasing with ^. This is 
expected due to the relation of the Polyakov loop and its 
conjugate to the free energy of a quark and an antiquark 
respectively [20j . In a system with more quarks than an- 
tiquarks it is relatively easy to screen a static antiquark 
by a quark, forming a virtual qq state, whereas a static 
quark can only be screened by a diquark, thereby forming 
a colorless three-quark state. 

At finite quark chemical potential, the order of the 
chiral phase transition can change. This is illustrated in 
Fig. [7] where the phase diagram of the PNJL model is 
shown. At low temperatures and large the transition 
is first order. The first order transition terminates at 
the tricritical point (TCP). With the actual value of the 
model parameters, summarized in Table. Illli the TCP 
appears at [Ttcp = 157, //tcp = 266) MeV. Beyond 
the TCP, at smaller chemical potentials, the transition is 
second order. This is consistent with the gross structure 
of the phase diagram, expected for QCD according to 
universality arguments [2l| . 

We explore the influence of the Polyakov loop on the 
thermodynamics, by considering observables that are re- 
lated to the conservation of the net-quark number such as 
the quark number density UqiT, fi) and the corresponding 
susceptibility Xq{T,fi). Both observables are obtained as 
derivatives of the thermodynamic potential f2 with re- 
spect to fi. The quark number density is obtained from 



9/i 



(3.1) 



With the thermodynamic potential of the PNJL model 
(|2?T5l) one finds 



6Nf 



(2^ 



9 



(+) 



(<I> + 2$e-^'"'/T^g-2B(+'/T 



-E^->/T 



$ + 2$e-^'"'/T_^g-2B<-)/T 



(3.2) 

The temperature dependence of the net-quark number 
density obtained with Eq. (|3.2p is shown in Figs. [5] and [H 
for various values of /i. The PNJL model results are also 
compared with those of the NJL model. Clearly there 
is a substantial change in rig when the Polyakov loop 
dynamics is included. 
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FIG. 5: Expectation values of the traced Polyakov loop $ (solid) and il> (dashed) in the chiral limit as a function of temperature 
T for finite quark chemical potentials. The lines in the left figure show $ and $ at ^ = 200 MeV and the chiral phase transition 
is of second order. The right figure is the results obtained at = 270 MeV, which correspond to the first order chiral phase 
transition. 
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FIG. 6: The dynamical quark mass M in the chiral limit as a 
function of temperature T for finite quark chemical potentials 
/i = 0, 200 MeV, 270 MeV from right to left. 
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FIG. 8: The net-quark number density riq/T'^ in the chi- 
ral limit as a function of temperature T at ^ = 100 MeV. 
The dash-dotted line denotes the result obtained in the NJL 
model. 
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FIG. 7: The phase diagram of the PNJL model in the chiral 
limit. The dot indicates the location of the tricritical point, 
located at {Ttcp = 157, /xtcf = 266) MeV. 



In the region above the chiral phase transition the stan- 
dard NJL model shows a relatively strong decrease of 
riq/T^ with increasing temperature. This behavior ex- 
actly reproduces the temperature dependence of a non- 
interacting gas of massless quarks, Uq = Nf{fj,^/Tr^ + 
/xT^). Thus, the leading term in Ug/T^ is proportional 
to fJ./T. Furthermore, the quark density increases as the 
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FIG. 9: The net-quark number density Uq/T^ in the chiral 
hmit as a function of temperature T. The lines, starting from 
the lowest one, correspond fi = 100, 150, 200 MeV. 



critical temperature Tc is approached from below. This 
is due to the decrease of the effective quark mass as the 
chiral symmetry is restored. The PNJL model shows 
a substantially different temperature dependence of the 
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quark density. First, there is a suppression of the one- 
and two-quark contributions to the density below T^, due 
to the interactions with the effective gluon field. Hence, 
the leading contribution to the net quark density is due 
to the three-quark states. Clearly this leads to a strong 
suppression of the quark density below Tc- Above Tc, 
the suppression is much less effective. However, because 
$ is still less than unity, a suppression compared to the 
quark density of a free gas at the same temperature and 
chemical potential remains. 

The quark number susceptibility measures the re- 
sponse of the quark number density to changes in the 
quark chemical potential. This observable is of partic- 
ular interest for exploring tricritical point. This is be- 
cause Xq is expected to diverge at TCP. This divergence 
is a remnant of the diverging fluctuations of the scalar- 
isoscalar sigma field [1^ [H, [2^ . The net-quark number 
susceptibility is defined by 



Xq 



9/i 



(3.3) 



In the PNJL model the dynamical quark mass M and 
the Polyakov loops $, i> implicitly depend on ^. Thus, 
besides an explicit ^-dependent contribution through Eq. 
p.3p there are also terms proportional to /^-derivatives of 
the effective condensates. In the PNJL model one finds 



,9$ 



- xf + xf ^ + xf^ + xf^ 

.(0) 



(3.4) 



where Xq corresponds to the 0-th order contribution 
given by 



6AV 
T 



(27r)3 



_ e(+) /t 

<I> + 4$e ^ 



3(l + l><I')e-2£<+>/T + 4<j,^ 



<i>e 



-4i5(+ 



$,$; 



(3.5) 



The functions A^'^) (M, $, T, /z) and the /x-dcrivatives 
of the condensates are introduced in Appendix |B] 

The quark number susceptibility at vanishing and at 
finite fi is shown function of the temperature in 

Figs. [To] and [U] for the PNJL and NJL models. At 
high temperature Xq approaches the ideal gas limit, 
Xq/T'^ — 2. In the chiral limit, the quark number suscep- 
tibility has a discontinuity at the chiral phase transition 
at finite u, as found in Landau- Ginzburg theory (see 
also [1^ [2y| for a detailed discussion of the quark number 
susceptibility in the NJL model). 

The effect of quark-gluon interactions on the quark 
number susceptibility, shown in Fig. IIO[ is similar to 
that discussed in the context of the net-quark density. 
There is a reduction of the quark fluctuations below Tc 
in the PNJL model relative to the fluctuations obtained 




0.25 0.5 0.75 1 1.25 1.5 



T/Tc 



FIG. 10: The quark number susceptibility Xq/T'^ in the chiral 
limit as a function of temperature T, at = 0. The dash- 
dotted line denotes the result obtained in the NJL model. 
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FIG. 11: The quark number susceptibility Xq/T'^ in the chiral 
limit as a function of temperature T. The lines correspond 
to the results in the PNJL model at ^ = 0, 100, 150, 200 MeV 
from below. 



in NJL calculations. The PNJL results for the T and 
/i dependence of Xq and Uq are in good agreement with 
recent LGT calculations of these quantities [13] ■ Thus, 
in contrast to the NJL model, the PNJL model provides 
a quantitative description of QCD thermodynamics near 
the phase transition. 

At the TCP the quark number susceptibility diverges 
as shown in the left panel of Fig. [121 We now analyze 
the different contributions in Eq. (|3.4|) . The right panel 
of Fig. [T2I clearly shows that the singular behavior of 
Xq is caused by the divergence of x)^^ ■ Both x^*-* and 
X^q^ have sharp peak structures at Tc, but, in contrast 

to ytq^'^ 1 they remain finite . 

In order to understand the relation between the singu- 
lar contributions to the susceptibilities, we consider an 
effective Lagrangian of the Landau type [2^ 

r!(T, /i) ~ l]o(r, + mM + ia(r, /i)A/2 + h){T, ii)M^ , 

(3.6) 

where a(T, /x) and &(T, \i) are temperature dependent co- 
efficients and TO is a source term, which breaks the sym- 
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FIG. 12: The quark number susceptibility Xq/T'^ the chiral hmit as a function of temperature T at the TCP 
The components xl^'^Xg*' ^-nd X?*"* a function of T at the TCP are defined in Eq. p.4p . 
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metry explicitly. At the end of the calculation we con- 
sider the "chiral limit" , and set m = 0. The equilibrium 
value of the order parameter M is determined by the gap 
equation dCl/dAI = 0. In the broken symmetry phase 
M = y^—a/b, while in the symmetric phase M — 0. 
The model exhibits a second order transition at a = 
for 6 > 0, and a tricritical point at a = & = [see e.g., 
Ref. [i^]]. The chiral susceptibility 



1 



dM 



dm a{T,n) + MPb{T,n)' 



(3.7) 



diverges along the second order transition line as well as 
at the TCP in the broken phase. 

On the other hand, the singular part of the quark num- 
ber susceptibility is given by 



— M — 



By using the gap equation one finds 



which implies 



Aq 



(3.8) 



(3.9) 



(3.10) 



where in the final step wc used the gap equation to 
eliminate M. Thus, the quark number susceptibility di- 
verges only at the TCP. The mixing of the chiral and 



quark number susceptibilities can also be interpreted as 
a consequence of ct-lu mixing at finite baryon density 



4. CHIRAL AND POLYAKOV LOOP 
SUSCEPTIBILITIES 

In the PNJL model the constituent quarks and the 
Polyakov loops are effective fields related with the order 
parameters for the chiral and Z(3) symmetry breaking. 
In LGT the susceptibilities associated with these fields 
show clear signals of the phase transitions. In this sec- 
tion we present the calculational scheme for computing 
the susceptibilities and discuss their relation to the spon- 
taneous breaking of the chiral and ^(3) symmetries. 

Consider first the generating functional of a scalar field 
theory 



Z[J] 



l?0exp 



i / d'^x {C[<j>] + J<j>) 



(4.1) 

where J is an external source. The second fimctional 
derivative of W[J] with respect to J yields the correlation 
function: 



S^W[J] 



,7=0 



(</.(x)0(2/))-(</>(x))(</>(y)). (4.2) 



SJ{x)SJiy) 

The effective action r[(f)] is introduced through the Leg- 
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endrc transformation 



r[0] = -W[J] - j d'^xJ{x)(j){x) . (4.3) 



The correlation function introduced in the Eq. (|4.2p can 
now be expressed by means of 



6''W[J] 



SJ{x)SJ{y) 



.7=0 



5(p{x)S(p{y) 



Thus, by using this identity, one can compute the suscep- 
tibiUty also from the effective potential, by taking deriva- 
tives with respect to the classical field 0. 

In the case of the PNJL model there are three different 
classical fields (f) — (M, $) that correspond to the or- 
der parameters; the dynamical quark mass, the Polyakov 
loop and its complex conjugate. Consequently, in order 
to compute the corresponding susceptibilities from the ef- 
fective potential, the relation (|4.4p must be generalized. 
We use the chain rule. 



6 SW 
5(f>j SJi 



SJk s^w 

56, SJkSJi 



i<j)j6(l)k SJkSJi ' 
(4.5) 

where J = {Jm ■, J<i>-, J^) is a vector composed of the 
source fields corresponding to the order parameters in 



It follows from Eq. (|4.5p that the susceptibilities are 
obtained by inverting the matrix composed of the second 
derivatives of the effective action with respect to classical 
fields. Following Ref. we introduce a dimensionless 
matrix 



C = 







^ml \ 




Clni 


Cu 


Cii , 


(4.6) 


C- 


C-n 


Cu / 





with the components 

1 d^^ 

1 d^n 



Cu 



1 d^n 



Cu = 

Cml 
Cml 



TA3 9$2 



Cu — Cii — 



rA3 9$2 ' 



Clm 
Clm 



1 d'^n 
1 d^n 



(4.7) 



TA2 dMd'P ■ 

Through Eq. (|4.5p a set of susceptibilities is defined by 



Xmm Xml Xml 
X = I Xlm Xll Xll 

/(.Im X.II 



(4- 



Here Xmm is the chiral susceptibility, as in the previous 
section, while xii Xu the diagonal Polyakov loop 
susceptibilities (see below). The off-diagonal terms cor- 
respond to mixed susceptibilities. 

In the pure gluon sector, the susceptibilities are related 
to the fluctuations of the $ and $ fields. Under the Z{3) 
transformation $$, $3, $3 and combinations thereof are 
invariant. Thus, the off-diagonal susceptibility Xu -^(3) 



(4.4) invariant, while the diagonal pieces Xu u ^^^t: 



Xu 
Xll 
Xu 



($$) _ ($)($) , 

($2) _ ($)2. 



(4.10) 



For vanishing quark chemical potential, xii ^-nd Xu co- 
incide but are not equal to XiT- In the low-temperature 
phase, where the Z{3) symmetry is realized, the Z(3) 
non-invariant components xii s-iid Xu strongly sup- 
pressed In the Z{3) symmetry broken phase at high 
temperatures, all components of the Polyakov loop sus- 
ceptibility are large. 

We also define the average susceptibility 



Xll = jixii +Xu + '^Xii)- 



(4.11) 



which corresponds to fluctuations of the real part of the 
Polyakov loop. This observable has been used in LGT 
calculations to identify the position of deconflnement 
transition in the QCD medium. 

Due to the presence of dynamical quarks in the PNJL 
model, the ^(3) symmetry is explicitly broken and the 
Polyakov loop is not a genuine order parameter. Never- 
theless, the fact that in LGT calculations the expecta- 
tion values of $ and $ are strongly suppressed in the low 
temperature phase [1^, suggests that the Z{3) symme- 
try is a useful guiding principle for constructing model 
Lagrangians (see also Figs. [3] and [5]). 



5. SUSCEPTIBILITIES AND THE PHASE 
TRANSITION IN THE PNJL MODEL 

Enhanced fluctuations are characteristic for phase 
transitions. Thus, the exploration of fluctuations is a 
promising tool for probing the phase structure of a sys- 
tem. The phase boundaries can be identified by the re- 
sponse of the fluctuations to changes in the thermody- 
namic parameters. In this section we focus on the phase 
structure of the PNJL model by studying the order pa- 
rameter susceptibilities, defined in the previous section. 



where Xij is given by the inverse of C, 



(4.9) 



In purely gluonic theory, where the Z(3) symmetry is exact, xil 
and Xu vanish identically in the low temperature phase. 
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FIG. 13: The chiral Xmm (dashed-line) and the Polyakov loop 
Xii (solid-line) susceptibilities in the chiral limit as functions 
of temperature T ior /j, = 0. 
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FIG. 14: The chiral Xmm (dashed) and Polyakov loop Xii 
(solid) susceptibilities in the chiral limit as functions of tem- 
perature T. The lines correspond to fx = (right), /j, = 200 
MeV (middle) and fj. = 270 MeV (left). 



5.1. Susceptibilities at vanishing quark chemical 
potential 



In Fig.[Tn]we show the chiral and Polyakov loop suscep- 
tibilities Xmm and Xil computed at /i = in the PNJL 
model in the chiral limit. The chiral susceptibility ex- 
hibits a very narrow divergent peak at the chiral criti- 
cal temperature Tch, while the Polyakov loop suscepti- 
bility shows a very different behavior: the peak is much 
broader and the susceptibility remains finite for all tem- 
peratures. The latter is due to the explicit breaking of 
the Z{3) symmetry by the presence of the quark fields 
in the PNJL Lagrangian. Nevertheless, Xii still exhibits 
a peak structure that can be considered as a signal for 
the deconfinement transition in this model. The peak 
position of Xu appears at T ~ 217 MeV below the chi- 
ral critical temperature Tch — 242 MeV. Thus, in this 
model the deconfinement phase transition occurs at a 
lower temperature than that of the chiral phase transi- 
tion. The separation between the two peaks is roughly 
25 MeV in the non-local model and about 8 MeV in the 
local model. 

Another interesting feature of Xii is the interference 
with the chiral susceptibility seen in Fig. [121 At the chi- 
ral transition, T = Tch, there is a narrow structure in 
Xii- We stress that this feature is not related with the 
deconfinement transition, but only due to the coupling to 
the chiral susceptibility. Thus, for the parameters used 
in the model, the deconfinement transition, signaled by 
the broad bump in Xih sets in earlier than the chiral 
transition at vanishing net quark density. 



As noted above, a similar procedure is usually applied in LGT 
studies of QCD thermodynamics. 



5.2. Susceptibilities at a finite quark chemical 
potential 

At finite chemical potential, there is a shift of the chi- 
ral phase transition to lower temperatures, as shown in 
Fig. [T] This is consistent with recent LGT results at 
finite quark chemical potential [sot . A lowering of the 
deconfinement temperature is also expected at non-zero 
net quark density. The position of deconfinement and 
chiral transitions can be determined by exploring the or- 
der parameter susceptibilities introduced in the previous 
section. With increasing chemical potential the temper- 
ature dependence of the Polyakov loop expectation value 
is flattening and for sufficiently large fiq it shows almost 
no variation with T. Consequently, the width of the 
Polyakov loop susceptibility is increasing with increasing 

In Figs. [U and [TS] we illustrate the temperature de- 
pendence of the susceptibilities for several values of the 
chemical potential. With increasing /i the peak position 
of the chiral and Polyakov loop susceptibilities are clearly 
shifted towards lower T and approach each other as seen 
in Fig. [141 At /zq — 185 MeV the two peaks coincide, 
which indicates that the chiral and deconfinement tran- 
sitions appear at the same temperature. For fi > fj,Q the 
Polyakov loop susceptibility Xii has a sharp peak at Tch 
and a broad bump above Tch. A similar behavior was also 
found in Ref. [3]. The peak at Tch in Xii is clearly due to 
an interference with the chiral phase transition, while the 
bump corresponds to the pseudo-critical point of the de- 
confinement transition. This structure is also seen in Xu 
at the tricritical point, located at (Tc = 157, Hc = 266) 
MeV. 

In the previous section we also introduced the aver- 
age susceptibility xu, which corresponds to fluctuations 
of the real part of the Polyakov loop. In Fig. [TSjthe tem- 
perature variation of xu for different values of the quark 
chemical potential is shown. It is clear from this figure 
that a behavior of xii is very similar to that of Xii- How- 
ever, the peak positions of the two susceptibilities differ 
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FIG. 15; The chiral Xmm (dashed) and Polyakov loop xu 
(soUd) susceptibihties in the chiral limit as functions of tem- 
perature T. The lines correspond to /i = (right), = 200 
MeV (middle) and = 270 MeV (left). 
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FIG. 17: The phase diagram of the PNJL model in the chiral 
limit. The solid (dashed) line denotes the chiral (deconfine- 
ment) phase transition respectively. The TCP (bold-point) is 
located at (Tc = 15, = 218) MeV. The results correspond 
to parameter set (b) in Table Hill 
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FIG. 16: The phase diagram of the PNJL model in the chiral 
limit. The solid (dashed) line denotes the chiral (deconfine- 
ment) phase transition respectively. The TCP (bold-point) is 
located at (Tc = 157, /Xc = 266) MeV. The parameter set (a) 
in Table Hill was used in the calculation. 



slightly. Thus, in the Polyakov loop sector the determi- 
nation of the transition temperature by using the sus- 
ceptibilities is not unique. Nevertheless, the differences 
are small, so in the following we use Xil to identify the 
deconfinement transition temperature in PNJL model. 

The peak positions of the Xmm and XiI susceptibilities 
are used to determine the phase boundaries in the (T, ji)- 
plane. We stress that the phase boundaries are strongly 
dependent on the model parameters. In Fig. [12] we show 
the phase diagram for the PNJL model obtained with the 
parameters from Table IH Clearly the boundary line re- 
lated with the chiral symmetry restoration, determined 
by the peak in the chiral susceptibility, coincides with 
that in Fig. [71 which was computed from the properties 
of the dynamical quark mass. Furthermore, the phase 
boundary corresponding to the deconfinement transition 
is identified with the position of the broad maximum ob- 
served in Xu- With these parameters the boundary lines 
of deconfinement and chiral symmetry restoration do not 
coincide. As anticipated in the discussion of Fig.[T31 there 
is only one common point in the phase diagram where the 



two transitions appear simultaneously. 

Recent LGT results both at vanishing and at finite 
quark chemical potential show that deconfinement and 
chiral symmetry restoration appears in QCD along the 
similar critical line 3l| . In general it is possible to choose 
the PNJL model parameters such that the critical tem- 
peratures of chiral and deconfinement transition coincide 
at /i = 0. The resulting phase diagram is shown in 
Fig. [T71 the parameters used in the calculations are sum- 
marized in Table UIH set (b). We note that this choice of 
the parameters is not unique. A shift in the position of 
the critical temperature at = can also be obtained by 
changing the To-parameter in the effective gluon poten- 
tial. Decreasing Tq from 270 MeV to 130 IVIeV results in 
Tch = 192 MeV, consistent with recent LGT calculations 

From Figs . [TBI and [T71 it is clear that in our model there 
is only a rather narrow region of finite where the decon- 
finement and chiral transition lines coincide. The slope 
of Tdcc as a function of /i is almost flat, indicating that 
at low temperature the chiral phase transition should ap- 
pear much earlier than deconfinement. So far there is no 
guidance available from first principle LGT calculations 
concerning the relation between deconfinement and chi- 
ral symmetry restoration at large values of the chemical 
potential. However, there are general arguments, that 
the deconfinement transition should precede restoration 
of chiral symmetry (see e.g. [H, HJ). In view of this, 
it seems unlikely that at T ~ the chiral symmetry sets 
in at the lower baryon density than deconfinement. In 
the PNJL model, the effective gluon potential param- 
eters were fixed by fitting quenched LGT calculations. 
Consequently, the parameters are taken as independent 
on ^. However, it is conceivable that the effect of dynam- 
ical quarks can modify the coefficients of this potential, 
thus resulting in /i-dependence of the parameters. Con- 
sequently, the slope of Tdcc as a function of ^ could be 
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FIG. 18: The diagonal xu ~ Xu susceptibility in the chiral 
limit as a function of temperature T for fi — 0. 
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FIG. 19: The off-diagonal Xmi = Xml susceptibility in the 
chiral limit as a function of temperature T for fi — 0. 



steeper Consequently, the cfFectivc Polyakov loop po- 
tential (|2.8p should, with ^-independent coefficients, be 
employed only for /i/T < 1. 

In Figs. [TH] and [11] we show further diagonal and 
off-diagonal components of the susceptibility x (|4.8[) at 
/i = 0. In this case the diagonal components xii ^-nd 
Xu coincide. These are the Z{3) non-invariant suscepti- 
bilities for <!> and $. These components are suppressed 
in the low-temperature phase, where the Z{3) symmetry 
is (approximately) restored. Around T ~ Tdoc the Xu,ii 
show a rapid drop, associated with the crossover transi- 
tion, where the expectation values and grow rapidly. 
In the high-temperature phase Xu u &re not necessarily 
suppressed since the Z(3) symmetry is explicitly broken 
there. The off-diagonal susceptibilities Xmi = Xmi ^^'^ 
also suppressed in the low-temperature phase and show 
a clear peak at Tch due to an interference with the chiral 
phase transition. These susceptibilities are non-invariant 
under both the Z{3) and chiral symmetries. Hence, they 
are suppressed well below the transition, where the Z(3) 
symmetry is restored and vanish above Tch, in the chi- 



Such a modification was explored in Ref. [sst where explicit fi- 
and Nf - dependence of Tq is extracted from the running coupling 
constant as, using the argument based on the renormaUzation 
group. 



rally symmetric phase. 

An alternative way to determine the dcconfinement 
and chiral transition temperatures, by locating the max- 
imum of the temperature derivative of the corresponding 
effective condensate, has been discussed in the literature 
(see e.g. ref. [l3l)- The temperature derivatives of the 
condensates can be expressed as combinations of the sus- 
ceptibilities with some T- and /i-dependent coefficients, 
as shown in Eq. (jB.3p . Consequently, in general the de- 
confincmcnt pseudo critical temperature obtained with 
this method does not agree with that obtained from the 
peak position of the corresponding susceptibilities. Only 
if the phase transitions are relatively sharp and the sus- 
ceptibilities show narrow structures can one expect that 
the transition temperatures determined using the differ- 
ent prescriptions coincide. In Fig. [^D] we show the deriva- 
tives of order the parameters at fx ~ for different values 
of To. 

In ref. Q the peak positions in the derivatives of the 
Polyakov loop and of the chiral order parameter coin- 
cide. Indeed, in the local version of the PNJL model with 
To = 270 MeV we reproduce this result, as illustrated in 
Fig. [20l However, for a slightly different value of Tq, 
the positions of the two peaks split. The parameters in 
the Polyakov loop sector were fixed from the lattice data 
in the heavy quark mass limit. This corresponds to a 
transition temperature To = 270 MeV in the pure gauge 
theory. In the presence of dynamical quarks the transi- 
tion temperature drops, thus it is not excluded that the 
value of Tq may depend on Nf, as noted in the introduc- 
tion. With To = 210 McV, the local PNJL model yields 
a splitting of the peaks in the derivatives of the Polyakov 
loop and the chiral condensate. Thus, the locations of 
the pseudo-critical points are strongly dependent on the 
model parameters. In the non-local PNJL model con- 
sidered in this work, there is no coincidence between the 
peak positions and the peak of |y is located at a lower 
temperature than that of ^y- for all values of Tq < 300 
MeV, as illustrated in Fig. [^Ok and c. We note that in 
the non-local PNJL model the chiral transition temper- 
atures obtained, on the one hand, from the derivative of 
the order parameter and, on the other hand, from the 
corresponding susceptibility are almost identical. Thus, 
at least in the chiral limit, the two methods for identi- 
fying the chiral transition, are equivalent. However, for 
the Polyakov loop, the difference is '--^ 10 MeV at fi ~ 0. 



5.3. Effective potential constraints 

As shown in Fig. I18[ xu is negative in a broad tempera- 
ture range above Tch. This is in disagreement with recent 
lattice results, where xii is always positive in the pres- 
ence of dynamical quarks (sij . Furthermore, the relation 
of the Polyakov loop susceptibility with the free energy 
of static quarks is also incompatible with negative values 
of Xu ■ A possible origin of this behavior could be the ap- 
proximation to the effective Polyakov loop potential used 
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FIG. 20: Derivatives of the dynamical quark mass dM/dT (dashed) and the Polyakov loop d^/dT (solid) in the chiral limit 
as functions of temperature T for = 0. 



in the Eq. 

The xii susceptibility can be related to the screening 
masses and m.i of the real and i mag inary part of the 
Polyakov loop correlation functions [36| : 



1 



1 



Xii 



In Ref. [36i the ratio 



!^^:.3 

rrir 



(5.1) 



(5.2) 



was found in the vicinity to the critical temperature. This 
result is qualitatively consistent with the perturbative 
calculation [13 ■ From the Eqs. (|5.ip and (|5.2p it is clear 
that the resulting xii should be positive at T = Tc. 

In Ref. [36| the screening masses and were calcu- 
lated with an effective potential that has the same poly- 
nomial form as used in the Eq. (|2.8p but with different 
values of the parameters. With the values given in table 
im effective Polyakov loop potential (|2.8[) yields 



(5.3) 



which implies a negative value of xu ^.t Tc- 

From this discussion it is clear that the behavior of xu 
depends crucially on the parameters used in the Polyakov 
loop potential. Thus, the constraints from Z(3) symme- 
try and from the lattice results for the equation of state, 
are not sufficient to warrant physically acceptable sus- 
ceptibilities. 



Recently an improved effective potential with 
temperature- dependent coefficients has been suggested 



m 



W($,$;r) a{T) 



b{T) In 



2 

1 - 6$$ 



4 ($3 + $3) -3($$)^ (^.4) 



where 



a{T) — ao + ai 



-02 



6(T) = 63 



(5.5) 

The polynomial in <& and <&, used in (|2.8p . is replaced by 
a logarithmic term, which accounts for the Haar measure 
in the group integral The parameters in (j5.4p were 
fixed to reproduce the lattice results for pure gauge QCD 
thermodynamics and for the behavior of the Polyakov 
loop. These parameters are summarized in the Table HVl 
In Fig. [21] we show the xii susceptibility calculated with 
the potential of Eq. (|5.4p . It is clear from this figure 
that the improved potential yields positive values for the 
Polyakov loop susceptibilities . In addition the peak po- 
sitions of the Xu ^"^d Xii susceptibilities almost coincide 
if the effective Polyakov loop potential is parameterized 
as in Eq. (|5.4[) . The phase diagram calculated with an 
improved potential, shown in Fig. 1221 is similar to that 
obtained with the previous choice of the Polyakov loop 
interactions. Fig. [111 This is also the case for the suscep- 
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FIG. 21: The diagonal xii = XU susceptibility in the chiral 
limit as a function of temperature T for ^ = 0. The effective 
Polyakov loop potential (|5.4p was used. 
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FIG. 22: The phase diagram of the PNJL model in the chiral 
limit. The solid (dashed) line denotes the chiral (deconfine- 
ment) phase transition respectively. The TCP (bold-point) is 
located at (Tc = 164, = 271) MeV. The parameter set (a) 
given in Table Hill was used in the calculations. 



gram and the order of the phase transition, using mean 
field dynamics for different values of the parameters. The 
properties of thermodynamic quantities related with the 
quark degrees of freedom, like the quark number density 
and susceptibility, were analyzed in the vicinity of the 
chiral and dcconfincmcnt transitions. 

We introduced susceptibilities related with the three 
different order parameters in this model, and analyzed 
their properties and their behavior near the phase tran- 
sitions. We have shown that there are as many as nine 
susceptibilities that can be used to identify the phase 
structure of the model. In particular, for the quark- 
antiquark and chiral density-density correlations we have 
discussed the interplay between the restoration of chiral 
symmetry and deconfincment. We argued that in the ac- 
tual formulation of the PNJL model a coincidence of the 
deconfincment and chiral symmetry restoration is acci- 
dental. 

We found that, within the mean field approximation 
and with the present form of an effective gluon potential 
the correlations of the Polyakov loops in the quark-quark 
channel show an unphysical behavior, being negative in 
a broad parameter range. This behavior was traced back 
to the parameterization of the Polyakov loop potential. 
We argue that the .Z(A^)-invariance of this potential and 
the fit to lattice thermodynamics in the pure gluon sec- 
tor is not sufficient to provide correct description of the 
Polyakov loop fluctuations in the presence of quarks in a 
medium. We note, however, that the improved potential 
of ref . [13 yields a positive, i.e. physical, xiu in quahta- 
tive agreement with the LGT results. 
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tibilities Xmm and Xii which have the same qualitative 
structure for both effective potentials. 

6. SUMMARY AND CONCLUSIONS 

We have explored the thermodynamical properties and 
the critical behavior of a system that exhibits an invari- 
ance under Z{2)) and chiral transformations. As a model 
we have use an extension of the NJL model for quarks 
with three colors and two flavors, which are coupled to 
an effective gluon field described by the Polyakov loop. 
In this model (the PNJL model) a non-local interaction 
instead of the point-like four-fermion couplings is em- 
ployed. The PNJL model exhibits two essential features 
of QCD: a spontaneous chiral symmetry breaking and 
a "confinement" like property. Furthermore, due to the 
symmetries of the Lagrangian, the PNJL model belongs 
to the same universality class as that expected for QCD. 
Thus, such a model can be considered as a testing ground 
for the critical phenomena related to the breaking of the 
global Z{i) and chiral symmetries. 

Within the PNJL model, we discussed the phase dia- 
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APPENDIX A: FUNCTIONAL INTEGRAL AND 
PARTITION FUNCTION OF THE MODEL 

The Polyakov loop L is a complex 3x3 matrix and 
can be diagonalized as in Eq. (|2.7p . Thus, the thermo- 
dynamic potential (|2.15p is a functional of complex vari- 
ables. However, as a physical observable the imaginary 
part of f2 must vanish. In this Appendix we show that 
this indeed happens. We follow the method that has 
been used in the context of a strong coupling QCD and 
the matrix model (ssl . [39j . 
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We start with the partition function, 



It is quite transparent that the pure gluon part is real. 
Thus, we focus only on the quark part Sq. Performing 
(A 1) functional integral over fermion fields one gets. 



with the action S being divided into three pieces: 

S\ii,^,L\ = dr J d?xC[i},4),L\ 

= Sg [L] + Sq[4>, V^, L] + 5i,t [V-, i'] ■ (A.2) 
The Polyakov loop effective potential U is included in Sg. 



where D denotes the Dirac operator 

D{^, ip')=i>-M- 7o(m + iA^) . 



(A.3) 



(A.4) 



Summing up the Matsubara frequencies, pq = iujn — i(2n + 1)ttT the partition function is obtained as 



\nZg = J VL{ifi,ip')Y^ J -^^Ti In D{ip, If' 



(27r)3 
d^p r, 



2Nf I VL{f,^') I A^[Tr,ln(l + L(v,,^')e-^^'^') + Tr^n [l + {v, ip')e-^^''') 



(A.5) 



taking the trace over color, flavor and Dirac variables in the above equation the imaginary part / of In Zq is proportional 
to 



/(Af,$,|.;r,Ai) = 



VLif,^') I -J^|tan-i 



(2n) 



3Im$ 



1 + 3Re$ (1 + e-E'+'/T) g-^'+'/T + e-3E(+} 



IT 



tan 



-3Im$ ( 1 - e-^'^'/T^ ,-b<-Vt 



1 + 3Rc$ (1 + e-^<-'/T) e-B(->/T + ^-ZEi-) jT 



(A.6) 



where ip and ip' dependence were suppressed and Rc$ = Re$ and Im<I> = — Im$ were used. Now let us replace 
the variables (p, with Lp —ip and tp' — > —ip'. The SU(3) Haar measure 'DL{tp,ip') is unchanged under these 
replacements while Im$ changes its sign. Therefore, the first and second terms are separately odd under the change 
of group variables and vanish after the integration. 

The thermal expectation values of complex $ and $ are evaluated as 



$) = ^ y" VL{if, (^')e'5'+-^'"' (Rezq ■ Re$ - Imz, ■ Im$) , 

(Rczq ■ Rc$ + Imz, ■ Im$) , 



(A.7) 



with 



Z 



(A., 



It is clear from Eq. (jA.6p that the imaginary part of the potential vanishes at /i = 0. Thus, the difference between 
($} and ($) comes only from the non-vanishing Imz^ at finite fi. This can be also seen in the matrix model for color 
SU(3) symmetry 
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APPENDIX B: DERIVATIVES OF EFFECTIVE CONDENSATES 



111 this appendix we summarize the derivatives of effective condensates with respect to n and T. Taking the 
/i-dcrivativcs in the coupled gap equations (|2.18p - (|2.20p one gets 



dM 


= 1 


dfj. 


A 


9$ 


T 




A2 


(9$ 


T 


9/i 


A2 



^A/ Xram ' ^ ^'I' Ami + A * 
2 Xml + -^A^ Xll + J^A^ Xll I , 



T 
A^ 



T 

T 



(B.l) 



where Xij a'l'e defined in Section 3] and the functions A^^-* are introduced as 



A 



M 



T3 J (27r)3 Ep 



a: 



6Af/ /■ 



T3 7 (27r)3 



_£;( + )/T 



1 - 3$e 



{2 + 3$e-^'"'/T_g-3B(-)/T| 



(B.2) 



The required temperature derivatives of order parameters are directly obtained from the gap equations as 



dM 


T 




A 




T 


df 


A2 




T 


df 


A2 



Xmm + ^ * ^ A * 



( ^Af Xm^ + ^^<I. Xzr+^^J, Xuj J 



T 



T 



where the functions A^"^^ are defined as 
,(T) _ 6Nf f cPp Mppip) 



A 



M 



T4 J (27r)3 



(B.3) 



A 



(T) 



Ta62|, , QNf f (fp 



2 dT J (27r)3 

-2£;<") /T 



i?(-){2 + 3|.e-^'"'/^-e-3^^'"'/T| 



i?(+){l-3$e-2^'*'/T_2e-3^'^'/Tj 



187V/ /■ d^p 



r3 y (27r)3 



5 



(+) 



5 



(B.4) 



APPENDIX C: PARAMETERS USED IN THE PNJL MODEL THERMODYNAMICS 



The compilation of parameters used in the model calculations is given in the following tables: 
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input 




= 92.4 MeV 


m-n = 135 MeV 


M 


= 335 MeV 


non-local NJL model 


A = 


684.2 MeV 


GsA^ = 4.66 


mo 


= 4.46 MeV 


local NJL model 


A = 


625.1 MeV 


GsA^ = 4.38 


mo 


= 5.31 MeV 



TABLE 1: Set of parameters for the NJL sector [l9|. The parameters of non-local NJL model were fixed for a = 10. 



Oo 


ai 


a2 


as 


63 


64 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



TABLE II: Set of parameters for the Polyakov-loop effective potential 0]. 
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set (a) 


set (b) 


A = 684.2 MeV 


A = 684.2 MeV 


GsA^ = 4.66 


GsA^ = 4.05 


To = 270 MeV 


To = 225 MeV 


rch(/i = 0) = 242 MeV 


rch(M = 0) = 180 MeV 


Tdccifi = 0) = 217 MeV 


Tdcc(Ai = 0) = 180 MeV 



TABLE III: Set of parameters in the chiral limit used in this work and tlie resultant phase transition temperatures, a was 
fixed to be a = 10. 



0,0 


ffli 


0.2 


h 


3.51 


-2.47 


15.22 


-1.75 



TABLE IV: Set of parameters for the improved Polyakov-loop effective potential [l^ . 
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